In the Claims 

1. (Original) A method for fast global motion estimation, the global motion 
defined by a plurality of parameters in which each parameter has an interval, the 
method comprising: 

a. providing a first and a second image, 

b. providing an initial estimate of each of two translation parameters, and 

c. determining the relative global motion between said first and second 
images using a symmetric gradient approach in an iterative process starting 
with said initial estimate of said two translation parameters, whereby said 
symmetric gradient approach provides the center of each parameter 
interval and results in improved global motion estimation convergence 
properties. 

2. (Original) The method of claim 1, wherein said step of providing an initial 
estimate of each of two translation parameters includes 

i. providing an initial interval for each said translation parameter, 
and 

ii. dividing each said translation parameter initial interval into two 
equal sub-intervals, 

and wherein said step of determining the relative global motion between said first and 
second images using a symmetric gradient approach includes 

i. providing a basic symmetric gradient formulation that includes 
said sub-intervals, 

ii. running in each iteration a point-wise linearization procedure 
on said basic symmetric gradient formulation, and 

iii. deriving in each iteration a symmetric linearization error based 
on said linearization procedure. 
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3. (Original) The method of claim 2, wherein said substep of providing a basic 

symmetric gradient formulation further includes using a motion parameters vector P 
representing the plurality of parameters. 



4. 



(Original) The method of claim 3, wherein said step of determining the 



relative global motion between said first and second images using a symmetric 
gradient approach further includes: for each iteration: 



second matrices, 

iii. calculating a vector Hflj 9 and 

iv. calculating a parameters vector p SGM using said combined 
matrix and said vector E?It, and using said symmetric linearization error for a 
continue/stop check. 



5. (Currently amended) The method of claim 4, wherein said substep of 



calculating respective first and second [H^Hj matrices includes calculating said 




calculating a combined matrix 




K Hz) using said first and 




matrices using respectively equations 59 [H* H 




dP k dPj 



and 




4 



6. (Currently amended) The method of claim 4, wherein said combined 
(HLK) matrix is calculated according to the following equation 62 

(g'<"4(M>(fl' fl );-) . 

7. (Original) The method of claim 1, wherein said global motion is selected from 
the group consisting from image translation, rotation, affine motion and panoramic 
motion. 

8. (Original) The method of claim 1, wherein said improved convergence 
properties include an improved convergence rate. 

9. (Original) The method of claim 1, wherein said improved convergence 
properties include an improved linearization error rate. 

10. (Original) A method for fast global motion estimation, the global motion 
defined by a plurality of parameters in which each parameter has an interval, the 
method comprising: 

a. providing a first and a second image, 

b. providing an initial estimate of each of two translation parameters, and 

c. determining the relative global motion between said first and second 
images using a bi-directional gradient approach in an iterative process 
starting with said initial estimate of said two translation parameters, 
whereby said bi-directional gradient approach provides the optimal 
location of each parameter interval and results in improved global motion 
estimation convergence properties. 
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11. (Original) The method of claim 10, wherein said step of providing an initial 
estimate of each of two translation parameters includes 

i. providing an initial interval for each said translation parameter, 
and 

ii. dividing each said translation parameter initial interval into two 
non-equal equal sub-intervals, 

and wherein said step of determining the relative global motion between said first and 
second images using a bi-directional gradient approach includes 

i. providing a basic bi-directional gradient formulation that 
includes said sub-intervals, 

ii. running in each iteration a point-wise linearization procedure 
on said basic bi-directional gradient formulation, and 

iii. deriving in each iteration a bi-directional linearization error 
based on said linearization procedure. 

12. (Original) The method of claim 11, wherein said substep of providing a basic 
bi-directional gradient formulation further includes using a motion parameters 
vector P_ representing the plurality of parameters. 



13. (Original) The method of claim 12, wherein said step of determining the 
relative global motion between said first and second images using a bi- 
directional gradient approach further includes: for each iteration: 

i. calculating separately for each of said first and second images 
respective first and second (MLK) matrices, 

ii. calculating a combined matrix ±£ BDGM using said first and 
second matrices, 

iii. calculating a vector H** DGM I t , and 
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iv. calculating a parameters vector P using said combined 
matrix and said vector H** DGM Ij t and using said symmetric linearization error 
for a continue/stop check. 



14. (Currently amended) The method of claim 13, wherein said substep of 
calculating respective first and second ^ H* h ) matrices includes calculating said 
matrices using respectively the following equations 78 and 79 

( H Hj =2^ 1 mid equation 

\==)kj dP k dPj 



15. (Currently amended) The method of claim 13, wherein said combined {I! 1 j£) 
matrix is calculated according to the equation &© h* dgm = h^ 1 H^ 2 . 



16. (Original) The method of claim 10, wherein said global motion is selected 
from the group consisting from image translation, rotation, affine motion and 
panoramic motion. 



17. (Original) The method of claim 10, wherein said improved convergence 
properties include an improved convergence rate. 

18. (Original) The method of claim 10, wherein said improved convergence 
properties include an improved linearization error rate. 
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